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ABSTRACT 

Changing the form of the star cluster mass function (CMF) can effectively change the 
upper end of the stellar initial mass function. The yields of from supernovae are very 
sensitive to the mass of the progenitor star. We show that by changing the parameters 
of the CMF, it is possible to change the yields of oxygen and magnesium by a factor 
of ~ 1.5 and of metals in general by a factor of ~ 1.8. 
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1 INTRODUCTION 

Virtually all stars form in clusters (Clarke, Bonnell & Hil- 
lenbrand 2000; Lada & Lada 2003). The type-II supernovae 
(SNII) of the most massive stars in these clusters are the 
main source of many heavy elements, in particular oxygen 
and a-elements. The yields of heavy elements from these 
SNII depend upon the mass of the progenitor (e.g. Tsu- 
jimoto et al. 1995; Woosley & Weaver 1995; Hoffman et 
al. 1999; Limongi, Straniero & Chieffi 2000; Rauscher et al. 
2002; Heger et al. 2002). Thus the yields of heavy elements 
will depend upon the initial mass function (IMF) of the most 
massive stars in clusters. Of particular interest in the study 
of SNII products are oxygen and magnesium that are both 
thought to be pure SNII products and which are relatively 
easy to observe. 

In large elliptical galaxies N/Fe, Na/Fe and Mg/Fe seem 
to be about a factor of 2 supersolar (Henry & Worthey 1999), 
as do Mg/Fe and O/Fe in stars of the Galactic halo and 
Thick disk (Prochaska et al. 2000). In the Galactic bulge, 
anomalies in ratios such as Mg/Ca suggest changes in the 
initial mass function (McWilliam & Rich 2004), which could 
also play a role in the a/Fe ratios, although their behaviour 
is usually attributed to differences in star formation his- 
tory. In this investigation we are mainly concerned with the 
overall oxygen and magnesium yields, which can be deduced 
from the peak in the abundance distribution function of stars 
(Pagel 1997); this peak for Mg is evidently at a somewhat 
greater abundance in the Galactic Bulge and the central re- 
gions of large elliptical galaxies than the approximately solar 
value that it has in the solar neighbourhood. 

In a recent paper Kroupa & Weidner (2003) showed that 
the (effective) upper IMF depends upon the cluster mass 
function (CMF). Low- mass clusters are unlikely to contain 
very massive stars as they do not have the gas reservoirs 
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from which to form them. A 1O 2 M0 cluster cannot contain a 
star > lOOM© and is unlikely to contain stars > 30 — AQMq . 
A 1O 5 M0 cluster, on the other hand, is very likely to contain 
some stars > lOOM© . This idea is supported by evidence for 
a truncated IMF in small clusters for stellar masses > lOAf© 
(Thilker et al. 2002). 

In the Solar neighbourhood the CMF is a power-law 
with slope ~ —2 between ~ 10 2 and 1O 4 M0 (Lada & Lada 
2003). There is evidence that the CMF varies with galaxy 
type (e.g. Kennicutt, Edgar & Hodge 1989; Thilker et al. 
2002; Alonso-Herrero et al. 2002; Youngblood & Hunter 
1999). 

In this paper we discuss the effect on the yields of oxy- 
gen and magnesium from SNII due to changing the CMF 
(and so effectively changing the upper-end IMF) . In Section 
2 we describe our method, in Section 3 we present the re- 
sults, and in Section 4 we discuss their implications for heavy 
element yields in galaxies of different types and masses. 



2 METHOD 

We select clusters from a cluster mass function (CMF) and 
then populate these clusters with stars from a standard IMF. 
Following Kroupa & Weidner (2003) we select star cluster 
masses at random from a power-law cluster mass function 
of the form N(M) oc M - ' 3 between lower and upper mass 
limits Mi ow and M up . 

Each cluster is then populated with stars drawn at ran- 
dom from a two-part Kroupa (2002) IMF of the form 



N(M) 



oc 



M -1 ' 3 0.08 < M/Mq < 0.5 
M -2 - 3 0.5 < M/Mq < 150 



such that the slope is Salpeter down to O.5M0 and then flat- 
tens. We ignore the brown dwarf regime as brown dwarfs, 
although numerous, contribute very little mass to the clus- 
ter. 

Stars are added to the cluster until the total mass of 
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stars exceeds the cluster mass. If the sampling produces a 
total cluster mass more than 2% greater than the desired 
cluster mass, then the cluster is completely repopulated. In 
each Monte Carlo run, the random sampling of clusters is 
continued until a total of 1O 9 M of stars has been formed. 

We record the numbers of stars greater than 8M and 
the enrichment from these stars. Initially very massive stars 
> 25Mq are expected to lose significant amounts of material 
through wind-driven mass-loss (e.g. Maeder 1992; Langer & 
Henkel 1995). For example, the yields quoted by Maeder 
(1992; table 6) for stars up to ~ 25M are very similar to 
those from Tsujimoto et al. (1995) and Thielemann et al. 
(1996). Above this mass, the yields drop sharply due to the 
lower mass of the eventual SNII. 

We take yields for oxygen (O) and metals (Z) from the 
combination of contents of stellar winds and final ejecta 
given by Maeder (1992) for solar-metallicity stars with 
strong mass loss (his Table 6). We assume that the yield 
of magnesium follows Tsujimoto et al. (1995) up to 25M , 
and that the Mg/O ratio does so above this. Between the 
tabulated points intermediate values are calculated by a lin- 
ear interpolation. In Fig. 0we plot the O, Z and Mg yields 
against the initial stellar mass. It should be noted that the 
assumption of these yields minimises the effects of the CMF, 
as it reduces the influence of very massive stars on the yields. 

Models of Galactic chemical evolution in the Solar 
neighbourhood have mostly either used a Salpeter IMF with 
an upper mass cutoff (e.g. Tsujimoto et al. 1995) or a steeper 
function such as the Scalo (1986) IMF (e.g. Chiappini, Mat- 
teucci & Gratton 1997). We select stars from a Salpeter IMF 
with an upper-mass cut-off of 15OM ; effective changes in 
the IMF are due solely to variations in the CMF. 



3 RESULTS 

In Table we show the effect of varying the parameters of 
the CMF on the number and mass of massive stars, and on 
the yields of O, Mg and Z. 

As reported by Kroupa & Weidner (2003), the slope of 
the high-mass end of the IMF depends upon the form of 
the CMF. The lower the mass of a cluster, the lower the 
probability of forming very massive stars. Indeed, low-mass 
clusters are unable to form very high-mass stars at all. This 
can be seen in the fractions of SNII progenitors with initial 
masses greater than 20, 40 and 70 M (F(> 20), F(> 40), 
and F(> 70) respectively). 

In Fig. [5] we compare the upper-mass IMFs of CMFs 
with (/3,M Iow ,M up ) = (2,50, 10 2 ) and (2, 10 3 ,10 6 ). For a 
CMF with parameters (2, 50, 10 2 ), which represents a CMF 
that is only able to form low-mass clusters, F(> 20) = 19%, 
F(> 40) = 3% and F(> 70) = 0.1%. For a maximum cluster 
mass of IOOM0 , it is almost impossible to form a star more 
massive than a few 10s of M . For a CMF with (2, 10 3 , 10 6 ), 
all of the clusters are massive. In this case, F(> 20) = 29%, 
F(> 40) = 10% and F(> 70) = 4%, which is close to what 
would be expected from pure random sampling from the 
Kroupa (2002) IMF without the constraint of fitting the 
CMF. 

The mass of stars formed per SNII changes between 
CMFs, rising from one per 1O8M for (2,50, 10 2 ), to one 



per 92M for (2, 10 3 , 10 6 ) — a rise of ~ 1.2 between the 
most extreme CMFs. 

It is a combination of the greater number of SNII per 
Mq of stars and the change in the fractions of SNII pro- 
genitor masses that produces a significant change in the 
yields between different CMFs. Between the (2, 50, 10 2 ) and 
(2, 10 3 , 10 B ) CMFs, the O, Mg and Z yields rise by a factor 
of 1.57, 1.40 and 1.83 respectively. 

Small differences in the yields can be seen if the slope of 
the CMF is varied. Between (1.5, 50, 10 4 ) and (2.5, 50, 10 4 ), 
the yields of O, Mg and Z change by a factor of 1.15 - 1.23. 

The total mass of stars forming was chosen to be 1O 9 M 
in order to allow a full sampling of the CMF. If the total 
mass of gas which forms stars is low, then different samplings 
of the CMF may have a small effect on the yields. For a CMF 
of (2,50, 10 4 ) sampled only from 1O 5 M of gas, the yields 
can vary by a factor of ~ 1.2. This variation comes solely 
from the random sampling of the CMF and the failure of 
some samples to include massive clusters and/or the success 
of others at doing so. 

Depending on the CMF, the number of SNII occur- 
ring per cluster changes by a huge factor, between 0.6 for 
(2,50, 10 2 ), to 75 for (2, 10 3 ,10 6 ). This would be expected 
to play a significant role in determining the effectiveness 
of the feedback of both energy and metals into the galaxy. 
When there are few SNII per cluster, it could be expected 
that feedback would be relatively inefficient. Much of the 
energy from the SNII would be used destroying the natal 
cloud, rather than used in spreading energy and metals far 
into the ISM. When the number of SNII in a cluster is high, 
it may be expected that the SNII would effectively form a 
superbubble, the influence of which may be widespread. 

The Solar abundance by mass of O is - 6x 10" 3 , and of 
Mg is ~ 7 x 10" 4 (Lodders 2003). The yields of oxygen are 
close to Solar, while the models we have used underestimate 
the yields of Mg compared to Solar. The exact details of 
the yields depend upon the models used. Whilst it would 
be expected that yields depend upon many factors such as 
the metallicity of the progenitor and the details of the mass 
loss, we do not think that the broad conclusions of this paper 
would change. Only in the highly unlikely case that yields 
scale with just the right power to offset the IMF could the 
yields be totally independent of the CMF. 



4 CONCLUSIONS 

We have shown that the yields of O and Mg may vary by 
factors of ~ 1.5, and of metals by ~ 1.8, depending on the 
form of the cluster mass function. This is due to the inability 
of low-mass clusters to form very massive stars (cf. Kroupa 
& Weidner 2003) thus changing the fractions of contributors 
of different masses to enrichment. 

In the Solar neighbourhood the CMF is a power-law 
with slope ~ — 2 between ~ 10 2 and 1O 4 M (Lada & Lada 
2003). There is evidence that the HII luminosity function 
(and so presumably the CMF) varies with galaxy type (e.g. 
Kennicutt et al. 1989; Youngblood & Turner 1999) with 
early-type spirals having fewer and less bright HII regions 
than late- type galaxies. This suggests a steeper CMF and/or 
a lower-mass cut-off (although this conclusion is debated by 
Thilker et al. 2002; their sample contains mainly late-type 




Figure 1. The yields of metals (dash-dot line), oxygen (full line) and magnesium (dashed line) against the initial stellar mass. Metal and 
oxygen yields are taken from Macdcr (1992) solar metallicity high mass-loss models, and the tabulated masses from Maeder are marked 
on the oxygen line. Magnesium yields are taken from Tsujimoto et al. (1995) and Thielemann et al. (1996) up to 25Mq, and are then 
taken to be 10% of the explosive yields of stars from Maeder (1992). 



Table 1. For a cluster mass function with slope (3 and lower and upper cut-offs M\ ovl and M up respectively, 
we show the mass of stars formed per SNII Mtot/WsNlli the average number of SNII per cluster Ns^n/N c i us , 
the fraction of SNII with masses > 20, 40 and 70M Q F(> 20), F(> 40), F(> 70), and the mass of oxygen, 
magnesium and metals produced per solar mass of stars formed Mq/M tot, A^Mg/^tot and Mz/Mtot- 



13 M low M up Mtot/WsNii N SNn /N clus F(> 20) F(> 40) F(> 70) M /M to t M Ms /M to t M z /M tot 
Mq Mq MQ cluster" 1 % % % xl0~ 3 xl0~ 3 xl0~ 3 
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galaxies, however). The upper limit of the HII luminosity 
function also varies with galaxy type, active star forming 
galaxies containing significantly more giant HII regions (e.g. 
Alonso-Herrero et al. 2002). 

It might be expected that dwarf galaxies have a CMF 
that has a low upper-mass cut-off. This is simply because 
dwarf galaxies have a far smaller reservoir of gas from which 
to form clusters, and the most massive clusters that they 
form would be expected to be fairly low. Thus dwarf galax- 
ies would be expected to have sub-Solar yields of heavy ele- 
ments. 

The Solar Neighbourhood has a CMF that is close to 
(2,50, 10 4 ) as shown by Lada & Lada (2003). This CMF 
produces close to the observed Solar abundances of O and 
metals, but under-produces Mg for these models. 

Massive elliptical galaxies are thought to form from ma- 
jor merger events (Sanders & Mirabel 1996). In such events, 
the star formation rate is thought to increase dramatically 
and super star clusters become a dominant mode of star 



formation. Thus, a CMF with high upper and lower-mass 
cut-offs may be a reasonable model for these galaxies. In 
such cases the yields could be a factor of > 1.5 higher than 
in dwarf galaxies, and > 1.1 higher than in normal spiral 
galaxies due solely to differences in the cluster mass func- 
tions. 
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Figure 2. A comparison of the upper-end IMFs produced by CMFs with (/3, M low , M up ) = (2, 50, 10 2 ) (lower curve) and (2, 10 3 , 10 6 ) 
(upper curve). The (2, 10 3 , 10 6 ) CMF clearly prod uces far more stars at high masses. 
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